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Abstract 

Gaseous Bose-Einstein condensates (BECs) have become an important test bed for 
studying the dynamics of quantized vortices. In this work we use two-photon Doppler 
sensitive Bragg scattering to study the rotation of sodium BECs. We analyze the 
microscopic flow field and present laboratory measurements of the coarse-grained 
velocity profile. Unlike time-of-flight imaging, Bragg scattering is sensitive to the 
direction of rotation and therefore to the phase of the condensate. In addition, we 
have non-destructively probed the vortex flow field using a sequence of two Bragg 
pulses. 
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1 Introduction 



Vortices are a hallmark of superfluids and have applications throughout the 
study of fluid mechanics and condensed matter physics [1,2,3]. Recently, gaseous 
Bose-Einstein condensates (BECs) under rotation have become an important 
test bed for predictions of the behavior of quantized vortices, which form 
highly regular lattices [4,5,6,7]. In these gases, time-of-flight (TOP) imaging 
has been an indispensable tool for observing the rotation of the cloud through 
the detection of the vortex cores. However, this technique only measures the 
superfluid density, and not the phase of the macroscopic wavefunction. Phase 
measurements would provide a more complete picture of vortex states, espe- 
cially useful in cases where individual vortices might not be easily detected. 
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In atomic gases, one can apply powerful optical and spectroscopic techniques 
[8,9,10,11,12] that afford new possibilities for performing phase measurements. 

In this work, we explore two photon Bragg scattering as a tool that provides 
information about the phase of a rotating Bose- Einstein condensate. We focus 

on both the microscopic and macroscopic signatures of the vortices in section 
2 of the paper, and present experimental data using this technique in section 
3. 



2 Vortex Lattices and the Bragg Method 

2. 1 Velocity Field of Rotating Condensates 

Vortex states in BECs have gained interest in recent years. Here we will give 
a brief introduction to their properties, and refer the reader to useful review 
articles by Fetter and Kevrekidis for further details [13,14]. Our starting point 
is a gas of bosonic atoms for which the order parameter ^(x, t) satisfies the 
Gross-Pitaevskii equation [15] 

where ^ = ^ is the reduced Planck constant, M is the atomic mass and 
Vext is the sum of all external potentials including the trapping potential. \E' 
is a single-particle Schrodinger wavefunction for which / l^pii^x = N, the 
total atom number. It approximates the many body physical description in 
the limit of weak interactions between the atoms and very low temperatures, 
both of which can be satisfied in the laboratory for a wide range of conditions. 
Since ^ is in general complex, it may be written as \Jn{x, t)e*'^('^'*\ where n 
is the superfluid density and S the phase of the wavefunction. The superfluid 
velocity field v = ^gV^ is proportional to the gradient of the phase. Since 
phase is only defined modulo 27r, it is readily apparent that the circulation 

§v ■ dl = m X must be quantized, where m — 1, 2, 3 The simplest state 

for which this condition is satisfied is a single vortex state with m — 1. 

For our purposes, we will assume that we have a harmonically trapped Bose- 
Einstein condensate in the oblate potential of a "TOP" , or time-orbiting po- 
tential magnetic trap. This trap has azimuthal symmetry in the x — y plane. 

It consists of a rapidly rotating magnetic bias field superimposed on a static 
quadrupole magnetic trap. The net effect is to move the magnetic field zero 
outside of the cloud, thus preventing nonadiabatic spin flips [16], and resulting 
in a time-averaged potential Vext — \M (cuf (x^ + y'^) + ujIz^). Within this po- 



2 



tential, a singly quantized vortex line aligned with the z direction has a phase 
= 0, where =tan~^(?//x) is the azimuthal angle. Thus the velocity field 

, ^ ^V7c ^ yx-xy 
M^^y) = M^S = --^ (2) 

where h is the Planck constant and M is the mass of the atom. 




Fig. 1. Microscopic velocity field of a vortex state. Velocity field of a single vortex 
is shown on the left. The arrows indicate the direction of the local velocity vector 
vi{x,y) whose magnitude is proportional to the arrow length. The projection vix 
as measured by Bragg scattering is plotted on the right. 

In Figure 1, we show the vector velocity field of a single vortex from Eqn. 2, as 
well as its x-projection vi^j. = vi-x as a 3-dimensional surface. From Eqn. 2 and 
the plots, one can see that although the vector magnitude It'll has azimuthal 
symmetry, the projection vix clearly does not - the x-velocity is largest along 
the y-axis and goes to zero along the x-axis. The region p = y/aP~+lp < ^ 
that contains the singularity at the origin has been removed from the plots 
for clarity. 

For large angular momentum, the superfluid typically breaks up into a lattice 
of singly-quantized vortices that form a triangular structure [5]. We have 
shown the corresponding column density profile of the atoms, proportional to 
/ l^pciz, in Figure 2a. In the vicinity of each vortex core, the particle density 
— s> within a region of radius ~ ^, where ^ is the "healing length", the mi- 
croscopic length scale that emerges from the GPE. Except at the condensate 
boundaries, the overall density distribution of the trapped gas is weh approxi- 
mated by the Thomas-Fermi profile n{p, z) = I^E'oP = ^~ ^'^ / R^z~ / ^tf- The 
Thomas-Fermi radii are Rz and Rtf along the z and p direction, respectively 
[15]. Typically, = 0.5pm, Rtf = 37/im and Rz = I3pm for our experimental 
parameters. To generate the column density profile, we multiply the Thomas- 
Fermi solution by a variational function f{x,y) for the vortex density. We 
used f{x,y) = l\f=i ^f^' where Si = \l {x - XiY + {y - ViY/^ and the A^^ 
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Fig. 2. Microscopic velocity field of a vortex lattice. Shown are (a) column density 
profile in the trap (approximate solution, left) and TOF images (experiment, center). 
On the right, the velocity projection v^^. is shown for a vortex lattice according to 
Eqn. 3. (b) Spectral density of Bragg scattered atoms for —2.7 it 0.3 mm/s (left), 
±0.3 mm/s (center) and 2.1 it 0.3 mm/s (right), showing microscopic as well as 
macroscopic flow. 

vortices are located at = {xi,yi), with i = 1,...N^ [17]. This allows the 
density to vanish at each singularity and gives us an approximate solution to 
the Gross-Pitaevskii equation which is valid so long as the vortex separation d 
greatly exceeds the core size ^. For our parameters d ~ Rtf/ \I N^/t: = 10/im 
~ 20 X ^. Adjacent to the vortex distribution, we have also shown an ex- 
perimental image of the column density profile which we have observed in a 
time- of- flight image of a rotating BEG, which confirms the lattice structure of 
the vortices. This method involves turning off the trap at time to and imaging 
the column density after a variable time of free flight, and the image roughly 
corresponds to the velocity distribution of the atoms just before to- 

The corresponding velocity field of a rapidly rotating BEG can be approxi- 
mated as the sum over all vortices. This is a valid approximation in the limit 
where the vortex separation greatly exceeds the core size: 



VN{x,y) 



h 

M 



i=l 



{x - XiY + {y- ViY 



(3) 



In Figure 2a (rightmost graph), we have plotted ■ x for a hexagonal lat- 
tice containing several vortices. One can clearly observe that the background 
velocity (i.e. away from the singularities) increases roughly linearly with the 
^/-coordinate. That is, when one coarse-grains over the velocity fields of the 
individual vortices, the resulting field is that of a rigid body rotation = Qy, 
where fl is the rate of rotation of the lattice [18]. There are also intriguing 
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signatures of the core regions in Figure 2b, which we will explore below. 



2.2 Bragg Scattering from Rotating Condensates 



2.2.1 Technique 

Time-of-flight images such as those we have shown in Figure 2a typically 
measure only the column density distribution of the atoms. That is, the signal 
in the image is proportional to and does not directly measure the phase 
S of the wavefunction ^ . One technique for probing S is two photon Bragg 
scattering [9,10], which is sensitive to the velocity of the atoms, which is oc VS. 
The full details of Bragg scattering process are beyond the scope of this paper, 
and therefore, we will concentrate only on a few key features which are relevant 
to our data. More details on Bragg scattering can be found elsewhere [19,20]. 

In brief, the Bragg method employs two laser beams with frequencies uj and 
uj + 2t:5. An atom scatters a photon from one laser beam into another. The 
net result is to impart to it a momentum q. For counterpropagating laser 
beams of wavelength \, q = |q| = 2/i/A is twice the momentum of a single 
photon. The two photon process and the energy-momentum relation are shown 
schematically in Figure 3a. For a condensate, q/M is typically much greater 
than the initial velocity of the atoms, and therefore, the diffracted cloud can be 
easily distinguished from the non-diffracted atoms. This is because the former 
have traveled an additional distance ^ q/M x ttof during the time-of-flight 
ttof after the trap has been shut off. The two clouds can then be separated in 
the images, as we show below. There is a resonance in the scattering of light 
when 



which expresses the conservation of momentum and energy. In the above equa- 
tion, the second term is the recoil energy, which must be provided by the en- 
ergy difference between the two photons. For sodium atoms near the principal 

2 _ 

resonance, 2^ = 100 kHz. The third term is simply the Doppler shift, which 
makes the Bragg technique velocity sensitive. It is this term which is of pri- 
mary importance to this work, as the Bragg process selects a group of atoms 
with the same projection of velocity Vx along the direction of the momentum 
transfer q — qx. For a trapped BEG, one also has to consider the effects of 
interactions. In the mean- field and local density approximations [20,21], this 



The amplitude of ^ does depend on the phase S, however, since the two are 
coupled during the time-of-flight evolution. 



5 



causes an extra 5mf = |^ frequency shift ^ , where ^ is the chemical potential 
This shift is a consequence of the Bogohubov dispersion relation for conden- 
sate excitations, as discussed in reference [9]. Free particles have 5mf = 0. 
In our case, for a stationary condensate, the Bragg resonance is peaked at a 
frequency 5o = + Smf ~ 101 kHz. Thus 6 = So + 



Energy 




(b) 



MamEntum 




Fig. 3. (a) Energy-momentum relation for atoms in a BEC showing the momentum 
transfer q from the two-photon process. The Bragg resonance occurs at a frequency 
difference 5 between the two laser beams that satisfies both energy and momentum 
conservation, (b) Diagram of experimental geometry. Vortices are created by rotat- 
ing a BEC about the z-axis by phase control of the transverse fields produced in 
coil pairs x and y that control the TOP (time-orbiting potential) trap. The Bragg 
beam containing frequencies iv and iv + 2715 is applied along the x-direction and 
retroreflected. 



2.2.2 Microscopic Flow Field 

Spatially selective Bragg scattering in principle allows one to "peer into" the 
complex velocity flow field of Eqn. 3. One may understand this in an intuitive 
way by considering the velocity surface Vx{x, y) in Figure 2. For a given Bragg 
beam detuning 5 we select a velocity through Eqn. 4. The region of space 
containing atoms that are moving at that velocity is given by the solution 
to the equation Vx{x,y) = Vx,o- Geometrically, this represents the intersection 
of the surface with a plane at height Vx,o, and is simply a velocity contour 
map. In Figure 2b, we show the resulting spatial distribution of atoms that 
are Bragg diffracted. The grey scale is proportional to the number of atoms 
within a velocity range of ±0.3mm/s, typical of the experimental resolution. 
Two features are clearly visible. For one, there is a broad horizontal band of 
atoms that are resonant with the laser whose center of mass y-coordinate yo 
varies as one varies the detuning of the laser beams, yo satisfies the rigid body 

^ The inhomogeneous density distribution also causes a broadening of the reso- 
nance. For our parameters, this is discussed in [22]. 
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condition Vx = ^y. We have varied the detuning to map out yo as a function 
of the velocity experimentally [22] . 

In addition, the spatial distributions contain information about the micro- 
scopic flow field of the vortices. The bands at constant yo in Figure 2b are 
not uniform; rather, there is a local structure to the velocity field associated 
with the fiow around individual vortices. This leads to beautiful patterns in 
the spatial profile of the diffracted atoms, including crescent and diamond 
shapes. These patterns are imprinted by the Bragg beams through the ve- 
locity selection. The existence of spatial signatures of the vortex state was 
first discovered by Blakic ct al. [23], who studied the Bragg scattering from a 
single vortex. They used a mean-field approach, incorporating the laser beam 
potential into the Gross-Pitaevskii equation to derive the spatial and temporal 
evolution of the Bragg scattered atoms. Observing signatures of these patterns 
in time-of-fiight imaging is complicated by the mean-field evolution during the 
time-of-fiight expansion, as discussed in [22]. In the following section we also 
discuss the necessary tradeoff between spatial and velocity resolution. 



2.2.3 Spectral Width 

The spectral width of the Bragg resonance is of critical importance, and is 
discussed in detail in a number of works. The finite temporal width tb of 
the Bragg pulse creates spectral broadening. Through the Doppler shift, this 
limits the velocity resolution Av^ to 

1 h 

= 5 

27rTB q 



A rotating BEG also has a spatial distribution of velocities. In order to achieve 
a spatial resolution Ax along the direction of the Bragg beams, the pulse 
duration must be short enough that the recoiling atoms from regions of the 
condensate that are separated by Ax do not spatially overlap during the pulse. 
Therefore, 

A. = ±rs (6) 



In our experiments we choose a time tb ~ 250/is. The recoil velocity is g/M = 
5.9 cm/s, while typical rotational velocities are Vrot ~ 3.5 mm/s and the 
Thomas- Fermi diameter is Dtp = 74/xm. This yields Av^ ~ Vrotf^O and Ax 
Dtf/5, which is a balance between good velocity and spatial resolution to 
observe the overall structure of the rotation. 

It is interesting to note the trade-off between high spatial and velocity resolu- 
tions, since their product AxAvx — h/M is a constant. For very short pulses 
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~ ^M/q = 8fis for our system, the vortex core regions could be resolved op- 
tically according to the above criterion; however, there will be a substantial 
spread of the velocities of the atoms which are Bragg diffracted. 



3 Experiment and Data 

In our experiments we produce a sodium BEC with typically 1 — 3 x 10^ 
atoms in a time-averaged orbiting potential (TOP) trap [16], according to 
the method described in [22]. Parameters for the trap are a radial gradient 
Bp — 12 Gauss/cm and a bias rotation of ujtop = 27r x 5 kHz. The measured 
transverse oscillation frequency is Up = 2% x 31 Hz, with Uz = \/8uJp. For 
details of experimental techniques, the reader is referred to several excellent 
review articles [24]. 

We produce the vortex-lattice by creating a rotating elliptical asymmetry in 
the horizontal x-y plane of the TOP trap [7]. The TOP trap employs a fast 
rotating bias field B = {Bx{t), By{t)) at a frequency ujtop much greater than 
that of the atomic motion. This creates a time-averaged harmonic potential. 
Therefore, we can create a slowly rotating elliptical potential by superimposing 
slow variations upon the fast oscillation of the bias field. To produce the 
fields, we combine the signals of two digital frequency synthesizers operating 
at frequencies uji = ujtop + ^ar and uj2 = ujtop ^^ar- These signals are each 
split, phase shifted and summed together to produce the two fields: B^it) ~ 
Bocos{u!it) + ecos{LJ2t) and By{t) = Bosm{u!it) — esin(a;2t), where e and cuar 
are the amplitude and frequency of the rotating asymmetry, respectively. The 
two currents are individually amplified using 100 Watt car audio amplifiers, 
and capacitively coupled to a pair of Helmholtz coils of approximately 10 cm 
diameter along the x and y directions, respectively (see Figure 3b). In order 
to maximize the number of vortices, we chose cuar = 27r x 22 Hz, which is 
very close to the frequency ~ O.Yujp that drives the quadrupole mode in our 
harmonic trap [25]. After applying the rotating asymmetry for 1.5 seconds, 
it was turned off and the atomic cloud allowed to equilibrate in the trap for 
another 1 to 1.5 seconds. This procedure reliably created vortex lattices with 
approximately 40 ± 10 vortices, as shown in Fig. 2a. 

3.1 Spatially Resolving the Velocity Field 

After producing vortices and allowing the lattice to equilibrate, we pulse the 
Bragg diffracting beams along the x-direction (see Figure 3)b) for a time Tb, 
while the atoms are still in the trap. The Bragg beams are detuned by 1.7 GHz 
from the F = 1 to F' = 2 resonance, and are created by back-reflecting a single 
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Fig. 4. Vortices probed by Bragg scattering. Outcoupled atoms (to the far right and 
left within each image) showed no particular structure for non-rotating clouds (a), 
whereas from vortices (b) and (c), the outcoupled atoms were tilted according to the 
direction of rotation. Images (a-c) were taken at 10 ms TOF. The tilt angle increases 
with respect to time of flight, as shown in (e). Each pair of Bragg frequencies is 
resonant with a thin strip of atoms parallel to the x-axis, as illustrated in (d). All 
images were taken at 5 = 102kHz. 

beam that contains two frequencies ujl and ul+^ttS. 6 is the difference between 
the frequencies of two rf synthesizers that are used to drive a single acousto- 
optic modulator. This creates 2 groups of diffracted atoms propagating to the 
left and to the right, respectively. We applied a Bragg pulse of square shape 
with tb = 250yus, and then turned off the magnetic trap within lOO/xs. The 
atoms expanded for a variable ttof before we took an absorption image using 
laser light resonant with the F = 1 — »• 2 transition in a 250yus pulse. The 
result is shown in Figure 4. We can clearly observe spatial structures in the 
outcoupled atom cloud arising from the rotation of the cloud. In figure 4a 
one can see the diffraction from an initially stationary condensate, and the 
diffracted (outcoupled) atoms appear to the right and left of the stationary 
condensate. No particular structure is visible. However, in Figure 4b, we have 
initially prepared a vortex lattice, which causes the diffracted atoms to form a 
tilted, elongated spatial pattern. Moreover, when we reversed the direction of 
the applied rotation (by replacing ljar —^^ar), the tilt angle with respect 
to the y-direction reverses, as shown in Figure 4c. 

We can understand our observations in terms of the coarse-grained velocity 
field discussed earlier: v = Q x f, with = proportional to the number 
of vortices. Since the Bragg process selects a group of atoms with the same 
Vx = ^y, the resonance condition is given hj 6 = 6o + 2yQ/X. Therefore, for a 
spectrally narrow Bragg pulse, with 6 < 6q and a counter-clockwise rotation, 
the resonance corresponds to a thin, horizontal band of atoms with y > for 
atoms which are Bragg scattered to the right, and ?/ < for atoms scattered to 
the left (the dark shaded regions in Figure 4d) . This band is identical to what 
is shown in Figure 4c, and within it there is a detailed microscopic structure 
near the vortex cores which is not resolved in our current experiment. As this 
band of atoms moves, the spread in velocities in the x — y plane causes part of 
the band to move up while another part moves down. Thus it forms a tilted 
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stripe whose angle increases with time, as observed in Figure 4e. At long TOF 
the stripe should become fully stretched along the vertical axis of the images. 



A calculation of the rotation frequency from the evolution of the tilt angle as 
function of time (see Fig. 4e), 6{t) = arctan(f2t), results in 1] = 27r x (15.4±1.1 
Hz). Moreover, one can use the locations of the diffracted cloud {Xfi,,Xl and 
Yr, Yl in Figure 4d) to spatially map the velocity field, which results in a 
similar value for the rotation rate [22]. This is in good agreement with the 
estimate based on the total quantized vorticity of the lattice. For that we note 
the fact that in the rigid body limit Q = {hNy)/{2m7rR'^) [18], and therefore, 
by measuring the number of vortices A^^, one can calculate the rotation rate. 
We used R = 37/im and was determined by a manual counting of the 
number of vortices from several images taken at long TOF. This resulted in 
Nv = 37± 7, which leads to 1] = 27r x (13.3 ± 2.6 Hz). 



Pulse 1 




Pulse 2 







V * V ' 


V1 V1-5V 



Fig. 5. In-Situ probe of the rotation. Two consecutive pulses separated by a hold 
time T can be used to non-destructively probe the rotating cloud. Shown are 10 ms 
TOF images after r = 7 ms (a), 5 ms (b), and 2.5 ms (c). The two diffracted groups 
appear as two tilted stripes on either side of the central, undiffracted cloud. Motion 
within the trap causes the first group of atoms (pulse 1) to change its tilt angle due 
to kinematic considerations, as explained in the text. 



3.2 Non-Destructive Probing of the Rotation 



The time-of-flight technique is a destructive method and cannot be used for 
in-situ observation of the dynamics of the rotating condensate. However, the 
Bragg method couples out a small fraction of the atoms, and therefore, can be 
used as a non-destructive probe. We demonstrate this by applying a sequence 
of two 0.5 ms duration Bragg pulses. Pulse 1 diffracted a group of atoms which 
then evolved for a variable time delay r in the harmonic trapping potential, 
after which a second group was created by pulse 2. The trap was then immedi- 
ately switched off, and an absorption image taken after 10 ms of time-of-flight. 
The results are shown in Figure 5 for 3 different hold times r. One can observe 
two diffracted stripes on either side of the undiffracted cloud in the center. 
These correspond to the two groups of atoms created by pulses 1 and 2. The 
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stripes tilted to the left are the result of pulse 2, and resemble those observed 
in Figure 4b. This indicates that the Bragg process is indeed non-destructive, 
and does not appear to significantly affect the rotation. Due to the confining 
potential, atoms are slowed down during the time r as they leave the cloud. 
Therefore, the diffracted cloud from pulse 1 has a lower velocity than that 
of pulse 2 when the trap is shut off. Therefore, it appears closer to the un- 
diffracted cloud after a long TOF than does pulse 2. Both pulses result in 
tilted images of the diffracted atoms; however, one can clearly see that pulse 
1 has a tilt angle and shape which is different from pulse 2. In fact, the tilt 
angle reverses sign after very short hold times r > 2.5 ms compared with the 
trap oscillation period of Tosc — 32ms. 

This rapid distortion of the diffracted cloud from pulse 1 within 2.5 ms evo- 
lution in the trap is somewhat counterintuitive. One might only expect the 
atoms to reverse their trajectory after Tosc/^ — 8 ms, that is, after a quarter 
period oscillation in the magnetic trap. This apparent paradox can be resolved 
when one considers the spatial extent of the diffracted cloud and simple kine- 
matics of the evolution. In Figure 5 we show the full potential V experienced 
by the Bragg scattered atoms. This consists of both the harmonic trapping 
potential and the factor of 2 larger mean-field repulsion from the undiffracted 
condensate at density tt-c, V = Vtrap + '^g^c [23]. Both atoms at the front end 
of the BEC as well as those at the back receive a kick and begin to move 
with a recoil velocity q/M which is considerably larger than the rotational 
velocity. Atoms at the front "chmb up" the potential for a time r = 2.5ms, 
and are therefore, slowed down, but atoms at the back end of the BEC do not 
change their energy substantially. The back end can thus catch up with the 
front end during the 10 ms TOF. Combined with the rotation of the strip, 
this causes a reversal of the tilt angle of the diffracted cloud from pulse 2. 
We have performed a numerical simulation of classical trajectories in a two- 
dimensional harmonic potential for our parameters. It confirms that the stripe 
should become vertical and tilt in the opposite direction for r > 2ms, which is 
consistent with our observations in Figure 5c. Therefore, the non-destructive 
method clearly has limitations, since the spatial TOF profile is the result of 
both the initial rotating distribution as well as its in-situ evolution during the 
hold time. 



4 Conclusion 

In conclusion, we have used Bragg scattering to directly measure the velocity 
profile of a rotating BEC. The technique is complementary to time-of-flight 
imaging, and has direct application to the study of non-equilibrium superfluid 
dynamics [1], where the tangle of vortices reduces their visibility. It might also 
be applied to measuring the normal fluid regime at finite temperature, where 
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vorticity is not present. Future experimental work will attempt to measure the 
microscopic flow fleld using Bragg diffraction. 

We thank Andrew Seltzman for experimental assistance. This work was sup- 
ported by the U.S. Dept. of Energy, the Army Research Office and by Georgia 
Tech. 
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